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ABSTRACT 



We study the collective kinklike normal modes of a system of several cylin- 
drical loops using the T-matrix theory. Loops that have similar kink frequencies 
oscillate collectively with a frequency which is slightly different from that of the 
individual kink mode. On the other hand, if the kink frequency of a loop is differ- 
ent from that of the others, it oscillates individually with its own frequency. Since 
the individual kink frequency depends on the loop density but not on its radius 
for typical 1 MK coronal loops, a coupling between kink oscillations of neigh- 
boring loops take place when they have similar densities. The relevance of these 



results in t he interpretation of th e oscillations studied by ISchrijver and Brown 



(I2OOOI ) and lVerwichte et al.l (l2004l ). in which transverse collective loop oscillations 
seem to be detected, is discussed. In the first case, two loops oscillating in an- 
tiphase are observed; interpreting this motion as a collective kink mode suggests 
that their densities are roughly equal. In the second case, there are almost three 
groups of tubes that oscillate with similar periods and therefore their dynamics 
can be collective, which again seems to indicate that the loops of each group 
share a similar density. All the other loops seem to oscillate individually and 
their densities can be different from the rest. 



Subject headings: Sun: corona-magnetohydrodynamics (MHD)-waves-scattering 



1. Introduction 



Transverse coronal loop oscillations w ere discovered by the Transition Region and Coro- 
nal E xplorer (TRACE) in 1998 (see, e.g. lAschwanden et all Il999l . I2OO2I : iNakariakov et al, 
I999I ). These oscillations were initiated shortly after a solar flare that disturbed the loops, 
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Since their first observation, transverse oscillations have been routinely observed an d stud- 



ied. Much before TRACE observations, the theory of loop oscillations was developed (ISpruit 



198ll : lEdwin and Robertslll983l : ICallylll986l ) and the different kinds of oscillations were stud- 
ied. The observed transverse motions have been interpreted in t erms of the fundamenta l 
kink mode of the fast magnetohydrodynamic (MHD) oscillation (INakariakov et al.lll999l ). 
which is the only mode that can produce the observed transverse loop displacement. 

In many cases the observed coronal loops belong to complex a ctive regions and are 



not is olated but forming bundles or arcades of loops. For example, in ISchrijver and Brown 



(|2000[) antiphas e tran sverse oscillations of adjacent loops were reported. In addition, in 
Verwichte et al.l (12004 ) phase and antiphase motions were observed in a post-flare arcade. 
On the other hand, it is cu rrently debated whether active region coronal loops are mo nolithic 



or multistranded (see, e.g. lAschwanden et al.ll2005l : iKlimchukI l2006l : iDeForestI 120071 ). In the 



multistranded model, it is suggested that loops are formed by several tens or hundreds of 
strands considered as miniloops for which the heatin g plasma properties are approximately 
uniform in the transverse direction (IKlimchukI l2006l ) . Most analytical studies about trans- 
verse loop oscillations have only considered the properties of individual loops. However, from 
the information provided by the observations, it is necessary to study not only individual 
loops but also how several tubes can oscillate as a whole, since their joint dynamics can 
be di fferent from that of a single lo op. Only a few works have considered composite struc- 



tures. 



Berton and Hevvaertd (119871) studied the MHD norma l modes of a periodic magnetic 



medium. iMurawskil (119931 ) and iMurawski and RobertsI (119941 ) studied numericall y the prop- 
agatio n of fast waves in two slabs unbounded in the longitudinal direction. In iDiaz et al. 



( 120051 ) the oscillations of the prominence thread structure were investigated. These authors 
found that in a system of equal flbrils the only non-leaky mode is the symmetric one, which 



means that all the flbrils oscillate in spatial phase with the same frequency. iLuna et al. 



( I2OO6I ) studied a system of two coronal slabs and found that the symmetric and antisymmet- 
ric modes can be trapped. A more complex system of two coronal cylinders was studied in 
Luna et al.l (120081 ). Four trapped normal modes were found and the interchange of energy 
between loops was shown by solving the time-dependent problem. In iTerradas et al.l (120081 ) 
a multistranded loop formed by ten strands was considered. The composite loop oscillates 
transversely as a whole with a global motion of the strands after an external disturbance. 
This work shows that the bundle of strands oscillates with a combination of collective modes. 
On the other h and, an analytical approximat ion to the normal modes of a loop pair has been 
carried out by IVan Doorsselaere et al.l (120081 ) . The authors assume the long wavelength ap- 
proximation and obtain an analytical dispersion relation for t wo d ifferent tubes together 
with the four kink mode polarizations described in iLuna et al.l (120081 ). 



In this work we aim to study the normal modes of a loop set with different physical and 
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geometrical properties by using the scattering theory. The scattering theo r y, or its matricial 



formulation called T-matrix theory (see, e.g. IWaterman and Truell 



1961 



Waterman 1969 



Ramml Il986l ). was first applied to magnetic tubes by iBogdan and Zweibell (Il985l ). These 
authors studied the interaction of acoustic plane waves with an ensemble of parallel magnetic 
fibrils distributed uniformly in the so-called spaghetti sunspot mode l. The authors derived 
and s olved the dispersion relation in the long wavelength limit. In iBogdan and Cattaneo 
( 119891 ) the frequency shifts and velocity eigenfunctions were calculated for the case of random 
fibril distributions of up to 100 flux tubes. Many other papers were published studying the 



cross section of a fib ril spot insonified by external acoustic waves (see lBogdan and Foxlll991 



Keppens et al.lll994j ). In all these papers a non- magnetized external medium was considered 
and the eigenfrequencies and eigenmodes of the acoustic oscillations were obtained. 

In this paper we generalize the method to a system with an external magnetized medium, 
in order to extend previous works to coronal loop conditions. Our model consists of an 
ensemble of parallel cylinders, without gravity and curvature. We consider uniform magnetic 
field in the internal loop medium and in the external or coronal medium. This assumption 
allows the interaction of the tubes through fast MHD waves. In addition, we explicitly 
calculate the eigenvalues and eigenfunctions of the normal modes of the model. 

This paper is organized as follows. In [|2]the loop ensemble model and the equations for 
its dynamics are presented. In ^ we briefly describe the T-matrix theory and apply it to 
our model. With this theory the exact eigenfrequencies and eigenmodes of two non-identical 
loops are investigated in §11 We study the dependence of the interaction with the relative 
density and radii of the loops. The study of three identical aligned loops is presented in ^ 
In the same section the interaction between three non-identical loops is considered. Finally 
in ®the results are summarized and the main conclusions are drawn. 



2. Theoretical model 

The equilibrium configuration used to model the loop set is a system of iV-cylindrical, 
parallel homogeneous flux tubes, with the 2;-axis pointing in the direction of the loop axes. 
All loops have the same length, L, and each individual loop, labeled j, is characterized by 
the position of its center in the xy-p\a.iae, rj = xje^. + yjGy, its radius, oj, and its density, 
Py The density of the coronal environment is po- The tubes and the external medium are 
permeated by a uniform magnetic field along the ^-direction (Bq = BqBz). The Alfven speed, 
va = Bq/ ^fjlp, takes the value fAj inside the j'-th loop and v^q in the surrounding corona 
(t;Aj < ^^Ao)- 
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Linear perturbations about this equilibrium for a perfectly conducting fluid can be 
readily described using the ideal MHD equations. In the zero-/? limit these equations can be 
written as 

^--iv^)pT = o, (1) 

where is the total pressure perturbation 

PT = —B,, (2) 

and Bz is the 2r-component of the magnetic field perturbation. The other perturbed quanti- 
ties, namely the velocity, v, the magnetic field perturbation, B, and the density perturbation, 
p, can be derived frompx- We have assumed a z-dependence of the perturbations of the form 
e~''''^^. The line-tying effect is incorporated by setting kz = qti/L where q is an integer num- 
ber. Hereafter we concentrate on the fundamental mode and take q = 1- We only consider 
problems for which the time dependence is a simple harmonic oscillation with frequency u. 
Then, the total pressure perturbation can be written in cylindrical coordinates as 

PT = e'^'^^--'^^P{r,^), (3) 

where iplr, ip) is a function that includes the dependence on r and ip. Inserting this expression 
in equation ([1]), we obtain the scalar Helmholtz equation 

Vlij + kl^l) = 0, (4) 

where _L refers to the direction perpendicular to the magnetic field Bq, i.e. to the 2;-axis, 
and kx_ is 

kl = ^^i^. (5) 
Hereafter, the _L symbol is dropped for the sake of simplicity. 



Normal modes 



The scattering theory, or its matricial formulation called T-matrix theory, provides an 
scheme to find analytically the norn aal modes of a system of scatterers in which waves are 
described by a Helmholtz equation (IRamml 119861 ). We fulfill the T-matrix theory require- 
ments because our ensemble of N loops can be considered a collection of scatterers and the 
perturbed total pressure is described by equation (jl]). 
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The T-matrix scheme states that the j-th flux tube generates an outgoing scattered 
wave, ipl^, in a field position r (in the two-dimensional xy-plane) that adds to the waves scat- 



tered from the other loops to produce the following net external field (iBogdan and Cattaneo 
1989h 



^(r) 



N 

E 

j 



(6) 



The scattered wave by the j'-th loop is produced as a response of an exciting wave produced 
by the external field minus the contribution of the mentioned loop, 



V'ix(r - rj) = '(/'(r) - ^J,(r - rj). 



(7) 



With equations and ([7]) the exciting field, ip^^, may be wr itten entirely in terms of th e 



scattered field, resulting in the self-consistency field equation (IBogdan and Cattaneo! Il989l ). 
This system of equations may then be closed by noting that the exciting and scattered fields 
are further related b y linear operators, T^, that describe the scattering proper ties of the 
individual fiux tubes (IBogdan and Cattaneo! Il989l : ! Wat erman! ! 1 9 69! : !Ramm!!l986! ) 



C(r-rj) = TVL(r-rj). (8) 

The key point is that the linear operators T-^ depend exclusively on the individual loop and 
external medium properties and can be directly computed through the boundary conditions 
on the loop-external medium interphase as we will see below. 

The external field to the j-th loop can be decomposed with equation ([7D as an excitation 
field on this loop and a scattered field by this loop. The excitation field has no sources in 
the j-th loop, i.e. it is the scattered field of the other loops, so it can be written as 



V'ex(^j,<^j) = Yl OimJra{koRj)e 



(9) 



where a^^ are the expansion coefficients of order m, that depend on the wave number 
and the frequency lj, and i?j and ipj are the local polar coordinates centered at rj, defined 
through Rj = |r — rj| and cosy^j = • (r — rj)/|r — rj|. Here Jm is the Bessel function of 
the first kind and order m and ko is k in the external medium calculated using equation ([5|). 
With equations ([8]) and Q, we find the scattered field in terms of an outgoing wave with 
sources at r 



J' 



(10) 



where T^^^ are the matrix diagonal elements of the operator T-i projected on the local basis, 
called T-matrix. The non-diagonal elements of this matrix are zero for axisymmetric tubes 



- 6 - 



( Bogdan and Zweibel 1985 ). The functions Hm^ are the Hankel functions of the first kind. 
With equations (Q, ([7]), ([9]), and (fTOj) . we find the following expression for the total field 



(11) 



in which the external field to 7-th loop is d ecomposed as an excitation on this loop and a 
scattered field by this loop ( jWatermanlll969l ). 



Following the development of iBogdan and Cattaneol ( 1l989l ). a linear algebraic system of 
equations for the complex coefficients al^ may then be obtained. We first substitute equation 
([6]) in equation ([7]) in order to obtain the self-consistency requirement 



(12) 



Next the exciting and scattered fields are replaced by their basis expansions, equat ions Qj 
and ffTOl) . and the translation formula (see appendix of iBogdan and Cattaneol Il989l ) is used 
to express the scattered wave centered in the i-th loop into an excitation at j-th flux tube. 
Finally, we obtain the following set of equations 



a- 



m + E E «n^nn^n-\n(ko|rj-ri|)e^(--)^J'=0, 
i^j n=-oo 



(13) 



where ipj is the angle formed by the center of the i-th loop with respect to the center of the 
j'-th flux tube. As we will see below with this equation we can find the al^ coefficients and 
the frequencies u, from which the spatial structure of the normal modes can be determined. 



From equation (jT3il . we see that the expansion coefficient of order m of the j-th loop, a^^, is 
coupled to all expansion coefficients of the other loops. This fact reflects the collective nature 
of the normal modes. With the al^ and equations ([6]) and (fTOj) we find the net external field. 

The internal or transmitted field is 



J2 PLJm{kjRi)e 



(14) 



where fcj is the transverse wavenumber inside the j-th loop calculated using equation ([5]). The 
Bessel functions of the second kind, Y^, are not considered in the expansion of the internal 



field (eq.[Tl]) because they are singular at the loop axes. The transmitted field (eq.[Tl]) can 
be calculated through the boundary conditions, namely the continuity of the total pressure 
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perturbation ([2]) and the radial component of the velocity at i?j = aj (see iGoedbloedlllQSSl ). 
In terms of the fields they are expressed as follows 

^Ufcj^j)|iJ,=a, = i^{koRi)\R,=a,, (15) 
K\kjRj)\R.=a. ^ V''(fcofij)|iij=aj . . 

where the prime is the derivative with respect to the function argument, iIj'{x) = dtp{x)/dx. 



Equation f[T^ is completely general for a system of cylindrical flux tubes (see lKeppens 



19951 . for non-axisymmetric expressions) and all the information of the individual loops is 
included in the T- matrix elements, T^^^. These elements are calculated through the bound- 
ary conditions at the interphase between the loop and the external medium. With equations 
( TTT]) . (IT^ . (ITBl) . and (ITBI) we find the T^^^ element expression 

_ fcj^fcoJm(fcjaj) J^m(fcoaj) - ko'^k^J'^{kfl^)J^{koaj) ^^^^ 
ko'^kjH^\koaj)J'm{k^aj) - k^'^koH''^\koa^)J^{k^a^) 



Equat ion (ITTl) is the generalization of iBogdan and Zweibell (119871 ): iBogdan and Cattaneo 



( 1l989l ) to the case of a magnetized environment. The zer oes of the d enominator correspond 



to the dispersion relation of the individual loop (see, e.g. ICallyl 119861 ). 



Finally, note that with the boundary conditions (eqs.[T5] and [TB]) it is possible to 
calculate the coefficients 

ni _ Jmikpaj) + Tl^H^\koaj) ■ 
which can be inserted into equation f|T4|) to obtain the internal field, ipl 



tr- 



Equation f[T^ is formally an infinite system of equations for an infinite number of un- 
knowns (ain). In order to solve it, we truncate the system into a finite number of equations 
and unknowns by setting ttmt+i ~ fo'^ azimuthal numbers greater than a truncation num- 
ber (m > rrit). To ensure the convergence of solutions, they must be independent of the 
truncation number rrit. With these considerations, the solution of equation (fT3|l reduces to 
solve a homogeneous linear system of N{2m^ + 1) equations and N{2mt + 1) unknowns. For 
this system of equations to have a non-trivial solution, its determinant must be zero. This 
requirement gives the dispersion relation as a transcendent equation. We solve the dispersion 
relation numerically and find the frequencies of the normal modes and with these frequen- 
cies we calculate the expansion coefficients. With equations ([6]) and f|T^ we find the net 
field in the external medium and the transmitted field in each loop. In all our calculations. 
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solutions are independent of the truncation number for values mt > 5 but we fix this number 
to rrit = 20 to more confidently ensure their convergence. With the method presented here 
we have obtained the results of the following sections. 



4. Interaction between two loops 

First, we compute the normal modes of two non-identical loops with the T-matrix 
theory outlined in ^ In this section we study the dependence of the interaction as a 
function of the density and radii of the loops. We consider a system of two loops with radii 
tti = a2 = a = 0.03L and separated a distance d = 3a. The first loop density is pi = 3po 
while p2 is changed from p2 = po to 5po to study its infiuence on the normal mode properties. 
We concentrate on the kinklike modes in which the individual loops move more or less as kink 
and suffer the largest transverse displacement. There are other higher order normal modes 
whose spatial structure is more complex, i.e. fiuting modes. We find four kinklike normal 
modes named P^, APy, Py, and AP^j., where P and AP refer to phase or antiphase motions 
of the loops, respectively, and the subscript x or y refers to the direction of the motion along 
the X- or y-axes. The frequencies of oscillation of these four modes as a function of P2/ po are 
displayed in Figure [H The bottom solid line is associated to the P^ and APy modes, which 
almost have the same frequency (see inbox Fig. [T]). The same behavior is found for the top 
solid line, which corresponds to the Py and AP^ modes. We see that the collective normal 
modes (solid fines) do not coincide with the kink frequencies of the individual loops (dashed 
lines), a discrepancy caused by the interaction between loops. This interaction is maximal 
when P2 = pi (dotted line), an d the normal mode s P^, APy, Py, and AP^ become the modes 



Sx, Ay, Sy, and Ay reported in iLuna et al.l (120081 ). The opposite situation takes place when 



P2 is sufficiently different from pi: the collective frequencies are closer to the individual 
kink frequencies and the system behaves as a pair of independently oscillating loops. In 
this regime, the P^ and APy modes correspond to the individual oscillations of the denser 
loop in the x- and ^/-direction respectively and possess identical frequencies, whereas the Py 
and AP^ modes are the individual oscillations of the rarer loop in the x- and ^/-direction, 
respectively, and also share the same frequency. Figure [1] can be interpreted globally as 
an avoided crossing of the kink modes of the loops: far from the coupling, each branch 
is associated to the individual loop kink mode, but near the avoided crossing motions are 
associated to the two loops to produce four collective kinklike modes. As for as kinklike 
solutions are concerned, loops interact for densities approximately in the range p2 = 2po to 
4po. 

The total pressure perturbation of the Px and AP^ modes is plotted in Figure [2] for two 
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Fig. 1. — Dimensionless frequency, ujL /vao, as a function of the internal density of the 
second loop. The bottom solid line is associated to the two kinklike normal modes and 
APy, which have very similar frequencies. In the same way, the upper solid line is associated 
to the Py and AP^ modes. In the inner plot a detailed view for p2 ~ pi is displayed showing 
that the solid lines are in fact double lines. The two dashed lines are the individual kink 
frequencies of each loop. We see that the kink frequency of loop 1 is constant and the 
frequency of loop 2 decreases with p2, because pi is constant and p2 changes. The vertical 
dotted line is plotted at P2 = Pi- Diamonds mark the frequencies of the modes represented 
in Fig. 



cases in which the loop interaction is important (p2 = 2.5po and 3.5po). The behavior of the 
other two modes, APy and Py, is analogous to that of the Px and AP^ modes and thus their 
spatial structure is not shown. In contrast to the case pi = P2, in which the in teraction is 



maxi mal and thus the two loops oscillate with equal amplitudes (see Fig. 2 of iLuna et al. 



20081 ) the solutions in Figure [2] display an imbalance in the oscillatory amplitude of the two 
loops. The largest amplitude of the pressure perturbation corresponds to the denser loop for 
the Px mode (see Figs. [2^ and[2)D), while it occurs in the rarer loop for the AP^ mode (see 
Figs. [2t and m). 

Secondly, we consider the same system of two loops but now the densities are fixed 
to pi = P2 = 3po, the radius of the left loop is ai = 0.03L, and the right loop radius, 
02, is allowed to vary. The distance between the tube centers is Sqm, where om is the 
averaged radius defined as om = (cti + 02) /2. With this condition the separation measured 
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(a) ojL/v^o= 2.189 (h) ojL/v^o= 2.054 




-4 -2 2 4 -4 -2 2 4 



-1.0 -0.5 0.0 0.5 1.0 

Fig. 2. — Total pressure perturbation of the fast collective normal modes Px and AP^ 
(plotted in the xy-plane) for a fixed density of the left loop (pi = 3po) and different densities 
of the right loop (p2)- The panels show the P^ mode with (a) p2 = 2.5po and (b) p2 = 3.5po; 
the APx mode with (c) p2 = 2.5po and (d) p2 = 3.5po. The frequencies of the modes 
are given on top of the corresponding panels. The dotted lines show the boundaries of the 
unperturbed tubes. Regions of the positive (negative) perturbed total pressure represent 
density enhancements (decrements), so that in (a) and (b) the loops move in phase in the 
x-direction, while in (c) and (d) they move in antiphase in the x-direction. 



in averaged radius units is constant. The frequencies of the four modes P^, APy^ Py, and 
AP^ are plotted in Figure [31 As in Figure [T] the collective frequencies (solid lines) are 
different from the individual kink frequencies (dashed lines), showing the collective nature 
of the oscillations. The chosen ra nge of radii are those me asured in TRACE observation 



of transverse oscillations (see, e.g. lAschwanden et al.ll2003l ). In Figure [3] we see that the 



collective frequencies are more or less constant; moreover the amplitude of the oscillation 
is more or less equal in each tube. Then, in the considered range of radii the interaction 



between kink modes is strong and does not significantly depend on the loop radii. 
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Fig. 3. — Dimensionless frequency, u)L/va.q, of tlie four collective kinklike modes Px, APy, Py, and 
APx (solid lines), as a function of the normalized right loop radius, 02/ L. As in Fig. [H the two 
individual kink frequencies are plotted (dashed lines), where the horizontal dashed line corresponds 
to the left loop and the one with a slight dependence on 02 corresponds to the right loop. 



5. Interaction between three loops 

5.1. Equal loop densities 

We first study the situation in which the density and radii of the loops is the same and 
find that there are eight kinklike normal modes, whose eigenfunctions are plotted in Figure HI 
with the modes ordered by increasing frequency. The lower frequency corresponds to a mode 
in which the three loops move in phase in the x-direction (Fig. H^), whereas in the higher 
frequency mode (Fig. Sh) the three loops move in phase in the |/-direction. This behavior is 
different from that of the system of two loops (see in which the higher frequency mode 
corresponds to the instead of the Sy mode. The modes of Figures H^, |Dd, H^, and|Di are 
kinklike while the other four modes of Figures Ht, HJi, H^, and Hf combine kink and fluting 
oscillations: the two left and right loops oscillate with a kinklike motion whereas the central 
loop oscillates with a fluting motion. We also refer to these modes as kinklike because at 
least one loop oscillates with a kinklike behavior. In these modes the central loop contributes 
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appreciably to the total field (eq.[TO]) with the multipole m = 2. Between the frequencies of 
the modes plotted in Figures Hli andH^ there are modes with the three loops oscillating with 
fluting motions and even with more complex structure associated to m > 2 solutions. Then 
we call these modes flutinglike. They are not further analysed because they do not produce 
transverse displacements of the loops. 

5.2. Different loop densities 

Now we consider the dependence of the interaction on the loop density. The loop radii 
are fixed to ai = 02 = as = a = 0.03L, the separation between adjacent loops is d = 3a, 
the densities of loops 1 and 2 are fixed to pi = 3po and p2 = 2po, and ps is changed from 
Po to 4po. Six kinklike normal modes, rather than eight, are found and their frequencies are 
plotted as a function of pa in Figure [51 There are six branches associated to the six kinklike 
modes, that have been labeled mi to uiq starting with the lowest frequency mode. We have 
chosen pi and p2 in such a way that loops 1 and 2 are basically decoupled (see Figure O 
is similar to Figure [1] and can be interpreted as two avoided crossings of the individual kink 
modes of the three loops. Far from the couplings, the loops behave independently. This fact 
is illustrated in Figure [61 In this figure we have plotted the modes associated to the branches 
mi, ma and ttiq in the top, central, and bottom rows, respectively. The modes m2, m4 and 
ms have an equivalent behavior, and have not been plotted. Far from the coupling region the 
mi and m2 solutions are associated to the individual kink oscillations of the denser loop in 
the X- and ?/-direction, respectively. In the same way, the branches ms and are associated 
to the individual kink mode of the intermediate density loop and the branches and mg to 
the individual kink oscillations of the rarest loop. On the other hand, at the couplings the 
loops interact by pairs as we see in Figure [3 The interacting pair oscillates with a collective 
normal mode whereas the other loop oscillates individually. In the first avoided crossing, for 
P3 = P2 = 2po, the branches ma and are coupled with ms and mg (see Fig. [5]), associated 
to loops 2 and 3, that oscillate collectively as we see in Figures [7t and [7^. The branches 
mi and m2 are uncoupled and loop 1 oscillates independently from the other two, as we see 
in Figure [7^. In the second avoided crossing at pa = Pi = 3po the branches mi and m2 
are coupled with ma and m4, while ms and mg are independent. Therefore, in this avoided 
crossing the interaction is between loops 1 and 3 (Figs. Uh and [711) and loop 2 oscillates 
independently (Fig. [7f). It is important to note that the collective modes of the two coupled 
tubes have a different frequency ordering with respect to the case of two loops, studied in 
§11 because of the presence of loop 1. 

Comparing Figures [H and [51 the coupling regions occur in a narrower range of density 
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Fig. 4. — Total pressure perturbation of the eight kinkhke collective normal modes of three 
identical loops. The densities are fixed to pi = P2 = Ps = 3po, the radii to CLl — 0,2 — O3 — 
a = 0.03L, and the separation between adjacent loops is d = 3a. 
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Fig. 5. — Same as Fig. [T]for the collective frequencies of three aligned loops plotted as a function 
of the density of loop 3. Solid lines correspond to the frequencies of the six collective kinklike 
modes. Dashed lines correspond to the individual kink frequencies of the loops. The horizontal 
bottom and upper dashed lines corresponds to the kink frequencies of loops 1 and 2 respectively. 
The other dashed curve corresponds to the kink frequency of the third loop, with variable density 
P3. Diamonds mark the frequencies of the modes represented in Figs. [6] and [71 

values in the three loop system than for two tubes. The physical meaning is that only loops 
with similar densities are coupled in the three loop ensemble. In this system, it is important 
to note that in the second avoided crossing at ps ~ pi loop 2 does not participate of the 
collective dynamics despite being the closest tube to the interacting loops. 

The results discussed so far in this subsection correspond to different densities of loops 
1 and 2. Nevertheless, if the densities of loops 1 and 2 are similar, their interaction is 
more important and the description of the dispersion diagram and the normal modes of 
the system is much more complex. In this case, there are eight kinklike normal modes. In 
this situation there are not only interactions between pairs of loops but also interactions 
between three loops. There are modes associated to the ensemble formed by tubes 1 and 2, 
individual oscillations of the cylinder 3 and the ensemble of the three loops depending on 
P3. A particular case of this situation is the three identical loops previously discussed ( §5.11) 
in which all modes are associated to the collectivity. 
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Fig. 6. — Same as Fig. Hlfor three values of ps far from the two couphng regions of Fig. [5l 
The densities of loops 1 and 2 are fixed to pi = 3po and p2 = 2po and ps is changed. The 
loop radii are fixed to = 0.03L. In the top, central, and bottom rows of panels 

the mi, ms, and mg modes are plotted for (a), (d), and (g) ps = 1.5po; (b), (e), and (h) 
P3 = 2.5po; (c), (f), and (i) ps = 3.5po. 



6. Discussion and conclusions 



In this work we have investigated the kinklike normal modes of a system of several loops 
with the help of the T-matrix theory. The results of this work can be summarized as follows: 



1. In the system of two non-identical loops, we have found four kinklike normal modes 
Px, APy, Py, and AP^. The frequencies of the Px and APy solutions are very sim- 
ilar as well as the frequencie s of the P„ and AP^ modes. This result agrees with 



Van Doorsselaere et al.l (120081 ) . who considered thin tubes (i.e. long wavelength ap- 



proximation). For fat loops the Px and APy modes, a s well the Py and APx, have 
different frequencies, as was shown in iLuna et al.l (120081 ). 



2. For a system of two loops we have investigated the dependence of the interaction be- 
tween kink oscillations as a function on the relative density of the loop pair. For 
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Fig. 7. — Same as Fig. Hlfor two values of ps at the maximal coupling (Fig. [5]). In the top, 
central, and bottom rows of panels the mi, m^,, and ttiq modes are plotted for (a), (c), and 
(e) p3 = 2po; (b), (d), and (f) ps = 3po. 

Pi = 3po we have found that the oscillations of the loops are coupled in the range of 
P2 between 2po to 4po and that the coupling is maximum at p2 = pi = 3po. Out- 
side this density range the loops are essentially decoupled and oscillate independently. 
This is qualitatively similar to the behavior of the anomalous modes described by 



Van Doorsselaere et al.l (120081 ). 



3. We have also studied the dependence of the interaction with the relative radii of the 
loops. We have seen that in the range of radii for which transverse loop oscillations 
have been observed the interaction depends very little on this parameter and the loops 
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strongly interact for all the radii considered. The explanation of this behavior is that 
in our loops the thin loop approximation can be applied and in this situation the kink 
frequency depends on the tube density and not on the radii. 

4. In the case of a system of three equal, aligned loops there are eight kinklike normal 
modes. The lower frequency mode corresponds to the three loops oscillating in phase 
in the x-direction, i.e. along the direction in which their axes are aligned, in agreement 
with the results of two identical loops. On the other hand, the upper frequency mode 
corresponds to the three loops oscillating in phase in the ?/-direction. This does not 
agree with the two identical loops situation, in which the upper mode corresponds to 
the two loops oscillating in antiphase in the x-direction. In fact, this property of the 
three-loop system is also true for ensembles of four or more aligned loops. 

5. We have made a parametric study of the kinklike modes in a system of three loops 
with equal radii and different densities by changing the density of loop 3, p^. We 
have chosen pi = 3po and p2 = 2po so that the interaction between loops 1 and 2 is 
negligible. We have found that the oscillations of loop 3 are coupled with loop 2 when 
Ps ~ P2, whereas loop 1 oscillates independently. Furthermore, loop 3 couples with 
loop 1 when p-^ ^ pi with loop 2 oscillating independently. If ps takes different values, 
the system is decoupled and the three loops oscillate independently. 



In this work, we have found that the interaction between loops regarding kinklike mo- 
tions depends strongly on the their individual kink frequencies. If these frequencies are 
similar, loop motions are coupled and the normal modes are collective. On the other hand, 
if the loop kink frequencies are quite different their motions are not coupled. Since the in- 
dividual frequencies depend on the loop density and radius, we have studied separately the 
influence of the two parameters. We have found that if the densities are quite similar, loops 
are coupled and the oscillations are collective. On the other hand, if the densities are quite 
different, the tubes oscillate independently. The range of densities for which the loops are 
coupled depends on the system properties and in the configuration of three loops this range 
is narrower than in the two tubes configuration. 

From the results s hown in t his paper we suggest t hat the antiphase motions reported in 
Schrijver and Brown I (120001 ) and lSchrijver et al.l (120021 ) are collective motions and, therefore, 
that the individual kink frequencies are similar but different from the collective observed 
frequency. If t he loop model p r esente d here is valid, both loop densities are also similar. 



In addition, in IVerwichte et al.l (120041 ) a loop arcade is studied and three groups of tubes 



oscillating with similar frequencies can be appreciated. The dynamics of each group of tubes 
can be interpreted as collective, although a detailed study of such configuration is needed to 
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relate the loop characteristics and the frequency of oscillation of the group. On the other 
hand, loops not belonging to these three groups do not share their frequencies with other 
loops, and so oscillate independently. This has to be interpreted as a sign that these loops 
have different d e nsitie s from those of the rest of the loops. It must be mentioned that in 



Verwichte et al.l (120041 ) all the oscillations are assumed as individual, but this is only true 
in the case of loops that do not share their frequency. If this assumption is applied to loop 
with a collective behavior it produces wrong results for the loop parameters. For example, 
if the loops actually oscillate with the lowest frequency collective mode, the assumption of 
these authors produces an underestimation of the magnetic field or an overestimation of the 
loop density. 

The T-matrix method shown in this paper can be easily applied to more complex con- 
figurations with gas pressure and tubes with flows or more complex systems of loops, i.e. 
arcades with myriads of loops or multistranded loops. It is expected that in such systems 
loops oscillate essentially independently except for loops with similar individual oscillation 
frequencies. 
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